Large Field Cutoffs Make Perturbative Series Converge 
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For A</>* problems, convergent perturbative series can be obtained by cutting off the large field configurations. 
The modified series converge to values exponentially close to the exact ones. For A larger than some critical value, 
the method outperforms Fade approximants and Borel summations. We discuss some aspects of the semi-classical 
methods used to calculate the modified Feynman rules and estimate the error associated with the procedure. We 
provide a simple numerical example where the procedure works despite the fact that the Borel sum has singularities 
on the positive real axis. 



1. INTRODUCTION 

It is intuitively clear that for Euclidean scalar 
field theory with quartic interactions, the sup- 
pression e "^^ makes the large field config- 

urations unimportant. This general idea has been 
tested explicitly by introducing sharp field cutoffs 
in models which can be solved numerically. A first 
example is the spectrum of the anharmonic oscil- 
lator which can be calculated very accurately [Q 
by using Sturm-Liouville theorem and requiring 
the wave function to vanish at some large value 
4'max rather than at infinity. A second example is 
the calculation Q of the zero-momentum correla- 
tion functions of Dyson's hierarchical model. This 
calculations requires the evaluation of the Fourier 
transform of the local measure and a sharp field 
cutoff can then be introduced at the first itera- 
tion, leaving the rest unchanged. In both exam- 
ples, if the field cutoff is taken large enough the 
errors due to the cutoff can be made exponentially 
small. 

In addition to being a natural procedure in nu- 
merical calculations, a large field cutoff cures the 
lack of convergence of the perturbative series. It 
removes the instabilities appearing at negative A. 
From a mathematical point of view, e"''"'^ has a 
Taylor expansion which is uniformly convergent 
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over the disk |(/)| < (f>max and one can then inter- 
change the sum and the integral. This statement 
has been checked numerically in Ref. |^ for the 
two cases mentioned above. It was also shown 
that the modifications to the Feynman rules due 
to the field cutoff could be calculated accurately 
by using semi-classical methods. In addition, the 
errors made at a given order, as a function of 
A and the field cutoff, can be estimated without 
knowing the numerical answer. At fixed A, this al- 
lows an optimal choice of the field cutoff. Beyond 
a certain value of A, this procedure outperforms 
Fade approximants and the Pade-Borel method. 
This is illustrated in Fig. |3| and discussed in Ref. 
1^. In the following, we first explain the gen- 
eral ideas with a simple example and then present 
some aspects not discussed in 

2. A SIMPLE EXAMPLE 

In order to understand the origin of the diver- 
gence and how to remedy it, we first consider the 
simple integral 
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If we expand e"''"'^ , the integrand for the order p 
contribution is e"'^^/^)'^ (p^P/pl and has its max- 
imum when 0^ — Ap. On the other hand, the 
truncation of e^'^'^ at order p is accurate pro- 
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vided that A^'* ^ p. Requiring that the peak of 
the integrand for the p-th order term is within 
the range of values of (f) for which the p-th or- 
der truncation provides an accurate approxima- 
tion, yields the condition A <C [l&p)~^ . One sees 
that the range of validity for A shrinks as one 
increases the order. We can avoid this problem 
by restricting the range of integration in Eq. (|l|) 
to |</>| < (jimax- As the order increases, the peak 
of the integrand moves across (pmax and the con- 
tribution is suppressed. Similar conclusions were 
reached in Ref. Q| using Lebesgue's dominated 
convergence. 

The coefficients of the modified series satisfy 
the bound \ap\ < -\/27r0^^^/p! and the modified 
series defines an entire function. However, we are 
now constructing a perturbative series for a prob- 
lem which is slightly different than the original 
one. This procedure is justified from the fact that 
we have an exponential control on the error SZ 
due to the restricted range of integration: 
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SZ < 2e-^*"- / d0e-(i/2)0^ (2) 

3. SEMI-CLASSICAL ESTIMATES 

Approximate bounds of the form of Eq. 
for lattice models or their continuum limit can be 
obtained by semi-classical methods. The general 
idea is that configurations with large field values 
are usually configurations of large action and they 
can be described by classical solutions with small 
fluctuations. A detailed calculation has been per- 
formed in Ref. ||^ in the case of the anharmonic 
oscillator. We used a dilute-gas approximation 
for configurations with one "lump" of large val- 
ues to estimate the effect of the field cutoff on the 
ground state. In the continuum limit, at imagi- 
nary time and at A = 0, the relevant classical 
solutions are 4'max^~^'^~'^°^ (in h = m = co — 1 
units). Such solutions can be observed in MC 
simulations for the harmonic oscillator. We first 
checked that proper thermalization occurred by 
comparing the distribution of the values of the ac- 
tion with a simple analytical result for quadratic 
actions as shown in Fig. [|. A typical configu- 
ration with large field values is shown in Fig. ^ 
together with the classical solution. 



Adapting standard semi-classical arguments, 
we obtained the zero-th order correction to the 
ground state: 

r+oo 

^ 4.-^/^0^,. / d^e-^' . (3) 

J 0max 

The effects of a non-zero A are obtained by re- 
placing 0"' by classical solution at A = 0. One 
obtains the approximate error at order n: 

\6Eo{X)\ ~ 5ij(°)e-(i/2)A0l„x + |a„+i|A"+i . (4) 

Both estimates are in very good agreement § 
with numerical results provided that (j)max > 2. 
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Figure 1. Action histogram for the harmonic os- 
cillator compared to {r[N/2])-^/^ S'^^^-^e"'^ for 
N = 400 sites. 
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Figure 2. A large field configuration (thin line) 
compared with 0maxe^'^^^°' (thick line) . 

The error due the field cutoff decreases when A 
increases, while the error due to the truncation of 
the perturbative series at a finite order decreases 
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when A decreases. One thus reaches a compro- 
mise at intermediate values of A where the total 
error is minimized. This is illustrated in Fig. |[ 
By increasing the field cutoff, the optimal value 
of A decreases. Using the approximate error for- 
mula, one can thus at the same time choose an 
optimal field cutoff and estimate the error. 
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Figure 3. Number of significant digits obtained 
with a large field cutoff (thick line) , for an expan- 
sion at order 3, compared to regular perturba- 
tion (thin line), Fade approximants for the regu- 
lar series (bottom set of dots), and the Fade-Borel 
method (top set of dots), at the same order. 

4. OTHER MODELS 

The generality of the basic principle invoked 
above (uniform convergence of the exponential 
on a compact neighborhood of the origin) im- 
plies that it can be used in many situations. In 
particular, for models where the Borcl sums have 
singularities on the positive real axis. A simple 
example is the integral M 
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^(1/2)0^ -|-((l+a)/3a)\/A</,^-(l/4a)A</,'' 
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Numerical results for a = 3/4 confirm that all 
the coefficients of the series are positive and that 
the Fade approximants of the series or of its Borel 
sums have all their poles on the positive real axis. 
The series defined with cutoffs in give results 
which are significantly better than regular pertur- 
bation theory when A > 0.1. This is illustrated in 
Fig. ^. Complementary discussions can be found 
inRcfs. 
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Figure 4. Number of significant digits for an 
expansion of order 3 of Eq. ^ for a = 3/4, 
with regular perturbation (thin line) and for a 
restricted range of (j> (thick line) . 



In lattice gauge theory, e.g. SU{2) in 4 dimen- 
sions, Ux,f_i — e'^3a.Af,.T .^j^jj contained in a 
sphere of radius — for each n . In the literature 
on lattice perturbation theory (with the excep- 
tion of van Baal), one usually replaces J dg^ by 
r^°° dA], since in the continuum limit the range 
becomes infinite. From our point of view, the 
lattice spacing a not only regulates the UV diver- 
gences but also the large order behavior of pertur- 
bation theory. We plan to treat this problem by 
following a path similar to scalar field theory but 
with different approximations to treat the mass- 
less quadratic theory with a field cutoff. 
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